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UTTOERWATER  EXPLOSION  BUBBLES 

I.   The  Effect  of  Compressibility  of  the  Water 
Joseph  B.  Keli.er  and  Ignace  I.  Kolodner 


1 .   Introduction 

As  soon  as  an  explosive  detonates,  it  is  converted  into 
a  gas  at  high  pressure.   In  an  underwater  explosion  this  gas 
is  called  the  "bubble".   It  expands  rapidly  until  its  pressure 
falls  below  that  of  the  suj:'roiinding  water  and  it  is  then 
compressed  again  to  a  high  pressure,  after  which  the  cycle  of 
expansions  and  contractions  is  repeated  periodically  with 
diminishing  amplitude.   Corresponding  to  the  pressure 
variations  in  the  bubble,  there  are  pressure  waves  transmitted 
outward  into  the  water.   The  first  pressure  wave  is 
a  shock.   These  waves   are    responsible  for  the  damage  caused 
by  the  explosion.   As  the  bubble  oscillates  it  also  rises  and 
flattens  due  to  gravity  and  the  influence  of  neighboring 
surfaces,  and  it  ultimately  breaks  up  or  vents. 

The  phenomena  described  above  have  been  the  subject  of 
extensive  theoretical  and  experimental  investigations  [1]. 
In  the  series  of  reports  of  which  this  is  the  first,  we  intend 
to  present  various  improvements  which  have  been  made  in  the 
theory  of  underwater  explosion  bubbles.   At  the  same  time  wo 
feel  that  these  reports  will  present  accurately  the  theory  of 
this  subject  in  its  present  form. 

The  simple  theory  of  underwater  explosion  bubbles  is 
based  on  three  main  assumptions:   1)  The  bubble  is  assumed  to 
be  spherical  at  all  times,  2)  the  pressure  is  assumed  to  be 
the  same  at  all  points  of  the  bubble  at  each  instant  of  time, 
3)  the  water  is  assumed  to  be  incompressible.   This  theory 
yields  two  somewhat  unsatisfactory  results:   1)  the  bubble  is 
foiond  to  perform  undamped  oscillations  of  constant  period, 
2)  the  pressure  variations  are  fo\ind  to  propagate  to  all 
points  of  the  water  instantly. 
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In  this  report  we  eliminate  these  features  by  modifying 
the  third  assumption  of  the  simple  theory  to  take  account  of 
compressibility  of  the  water.   We  then  find  damped  oscillations 
of  diminishing  period,  a'ld  finite  propagation  speed  for  the 
pressure  waves  in  the  water.   The  theory  and  results  are 
essentially  the  same  as  those  of  Zoller  [3],  with  some 
additions . 

In  section  2  the  problem  is  formulated  mathematically, 
and  this  formulation  involves  one  novel  feature  which  may  be 
of  more  general  interest.   This  is  the  use  of  the  wave 
equation,  rather  than  the  ii'.sual  Laplace  equation,  for  the 
potential  function  in  the  v/ater.   It  is  surprising  that 
although  the  wave  equation  is  closer  to  the  exact  non-linear 
equations  for  the  potential  than  is  the  Laplace  equation, 
still  it  is  widely  held  that  the  wave  equation  is  applicable 
only  to  small  motions,  iirhile  Laplace's  equation  is  employed 
for  large  motions. 

A  demonstration  of  the  applicability  of  the  wave 
equation  to  large  motions  has  been  given  by  G.  I.  Taylor  [2], 
who  compared  the  exact  solution  for  the  flow  produced  by  a 
liniformly  expanding  sphere,  with  the  approximate  solution 
given  by  the  wave  equation.   The  agreement  was  excellent, 
even  for  very  large  velocities  of  the  sphere,  while  the 
corresponding  solution  based  on  Laplace's  equation  misses 
some  of  the  main  feati.\rcs  of  the  flow.   3ccause  of  the 
similarity  of  Taylor's  problem  to  the  present  one  it  seeiis 
clear  that  the  wave  equation  should  also  yield  a  good  solution 
in  this  case,  and  thir.  cxi^ectation  is  borne  out  by  the  results. 

In  section  3  the  iiathematical  probl<.m  is  solved  and  in 
section  [|.  the  direction  field  which  characterizes  the  result, 
is  analyzed.   Sketches  of  this  field,  and  the  corresponding 
one  obtained  by  siooclalizing  the  result  to  an  incompressible 
fluid,  are  given.   A  graph  of  bubble  radius  versus  time  based 
on  a  numerical  integration  for  the  compressible  case,  is  also 
given  and  compared  with  an  experimental  curve.   Graphs  are 
also  pivcn  from  which  the  pressure  in  the  v;atcr  can  be 
determined  as  a  function  of  position  and  time. 
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In  section  5  an  analytic  solution  of  the  linearized 
equation  for  the  bubble  radius  is  given.   This  applies  to  a 
bubble  whose  pressure  is  not  much  different  from  the  surrounding 
hydrostatic  pressure, 

flection  6,  the  conclusion,  discusses  the  results  in  the 
light  of  other  work  on  this  subiect.   In  the  Appendi.^  the  wave 
equation  is  derived. 

2.   Formulation 

V/e  consider  a  sphere  of  gas  (the  bubble)  of  initial 
radius  A   surrounded  by  an  unbounded  fluid  (the  water) 
initially  at  rest.   VJe  assvime  that  the  bubble  remains  spherical 
at  all  later  times,  that  the  pressure  throughxout  the  bubble  is 
constant  at  each  Instant  of  time  It  and  that  the  pressure  TT(x) 
in  the  bubble  is  related  to  its  radius  A{t)    by  the  adiabatic 
relation 

(1)  TT(t:)  =  K(^  A^-r))'^ 

We  further  assume  that  the  water  velocity  is  derivable 
from  a  potential  function  £(R,tJ)  which  depends  only  upon  radial 
distance  R  from  the  bubble  center  and  time  Z  ,    and  satisfies 
the  wave  equation  ( see  Appendix) 

(2)  Al-^i^^=o   . 

Here  C  is  the  soimd  speed  in  the  fluid,  assumed  to  be  constant. 
The  pressure  P(r,T;)  in  the  water  is  given  by  the  Bernoulli 
equation 

(3)  P(r,tr)  =  Fq  -  p(i^+  ^I^)   . 

Here  p  is  the  density  of  the  fluid  .;nd  P^  is  its  initial 
pressure . 

At  the  bubble  surface,  pressure  must  bo  continuous  and 

the  rate  of  change  of  bubble  radius  must  equal  the  particle 

velocity  in  the  water.   Thus  \io   have 


ik)         K(^  A^(t))-Y.  p^_p(J^+  1  J;)    at   R=A(r) 

(5)  i  =  A^  at   R  =  A(  -r)    . 

The  problem  is  to  solve  (2),  (1^.),  (5)  for  '^{'R,r)     and 
A(  Z)     given  that 

(6)  A^(0)  =  Aq,  A(0)  =  A^;  J(R,0)  =  0,  R  >  A^;  i^(R,0)  -  0,  R>  A^ . 

In  the  simple  theory  C  is  infinite  so  that  -ij-  _^   in  (2)  does 

C    ^^ 
not  appear--the  problem  is  otherwise  the  same,   '/e  will  obtain 

the  solution  in  that  case  as  a  specialization  of  our  result. 

Before  solving  this  problem,  we  will  introduce 

dimensionless  variables  as  follows. 


where 


(7)   A=aL,  r=  Tt,   J=L^t"-'-<I>,   R=Lr,   C=LT"-'-c,   A^  =  a^L 


inev^)'''.  '   .    T  =  lJ;^ 


L  =  ,^r°vvr 


a  =a^, 


'•,'ith  these  changes,  (2),  (Ij.),  (5).  (6)  become,  v;ith 
> 

(8)  ^  ct)  -  ^  (j,^^  =  0 

c 

( 9 )  't'^  +  i-  a^  +  a"^^  =  1    at  r  =  a 
(IG)  '^r  ~  ^   at  r  =  a 

(11)  a(0)=a^,  a('c]=ai^i    C>(r,0)=G,  r>a^;  -j)^  (  r,  l)  =  0,  r  >  a^  . 

3 •   Solution 

To  solve  (3)-(ll)  we  choose  a  particular  solution  of  (£), 
namely 

(12)  ^-_£mt.rll     ,    ;;=  t-  r>   . 
Then  (7),  (10)  become 
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(13)      ci'^r'(2:(t,a)) +|a^+  a~3Y   ^  1 

ilk)  a  =  -a"^f(^(t,a))-(ac)''^f'(?(t,a)) 

Elliuinating  f'(^(t,a))    bctv;oon   (13),    {Ik)    and   differentiating 
the  result  with  respect   to  t  yields 


(15)  (a-c)f'(^(t,a))   i-Ui  a2+a-^Y.i_c^)a2]. 


0 


Nov;  clL-ninating  f'(^(t,a))    between   (13),    (1^)    and   sir.iiDlifying 
yields 

(16)  (a-  c)(+aa  +  ^  a'^ -a^^Y^.  1)  .  (j^2^  ^^^_  2)a-3Y+2)d  =   0      . 

Equation  (16)  is  a  non  linear  second  order  ordinary 

differential  equation  for  a(t),  and  the  initial  data  arc  given 

in  (11).   Once  a(t)  Ig  laiovm,  f  and  f  can  be  determined  from 
(13),  (lU),  which  yield 

P    2   .2 

(17)  f(^(t,a))  =-a^a  +  ^(^  +  a-^Y_i) 

f'(^(t,a))  =  -a(^+a-^Y  .3) 

Since,  from  (11),  ^{0.    ,0)    =  0  v.'c  must  have  f(0)  =  0,   Using 
this  in  (17)  yields 

(IQ)  a^  =  c  1-  Jc^  +  2(1  -a;3Y)    . 

As  a  mu3t  bo  less  than  c,  v;o  must  choose  the  minus  sign  in  (I8). 

TI^cn  a  ^  0  according  as  a  5  1.   There  is  no  real  value  of  a^ 
O  ^  <-'      o  -•  o 

/c^   A-1/3Y   •  • 
if  a  <  I -^  + 1  j     ,  i...  if  a   is  too  small,  in  v/hich  case 

the  present  theory  fails  since  the-  initial  bubble  pressure  is 
too  higj.1.   If  c  is  large  compared  to  2(1 -a'-'*'^)  then  (I8)  becomes 

-1-f  a-^Y 
(19)  i  ft^ r-^-    . 


l/c  note  that  f  becomes  -a  a  and  a  becomes  0  when  c  is  infinite. 
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Since  t  docs  not  occur  explicitly  in  (16),  V7C  may  reduce 
the  equation  to  first  oi'dcr  by  defining 

(20)  v(a)  =  a 

Then  (16)  becomes 

v.(a)  =  JlpLi^   , 
^     av(c-v)    * 

(21) 

h(v,a)  =  ^(v^  -3v^c  -2v  -2c)-a"^Y(3Y  _i)(v  -j^) 

In  terms  of  v(c.),  the  solution  is  given  by 

^a 
o 

Thus  wc  must  solve  (21)  for  v(a),  then  obtain  a(t)  from 
(22)  and  finally  determine  f  from  (17)» 

To  calculate  the  pressure  P(r,t)  v;c  have  from  (3).  (12) 

(23)         P(r,t)  =pji-{i;H.  i:^.££;. -£^)1   . 

2c  r   cr-^   2r^ 

The  argument  of  f  in  (23)  is  ^(t,s).   Hov/evcr  in  (1?)  f  is 
given  in  terms  of  argiuncnt  ^(t,a(t)).   To  find  f(^(t,r)),  wc 
define  f(^(t,a))  =  g(t),  and  g  is  then  given  by  the  right  side 
of  (17).   Then  if  we  choose  x  such  that 

(2li.)  a(:c)  -ex  =  r  -  ct 

wc  havo 

f(^,(t,r))  =  g(x) 

IIov;   (23)   roay  be  rewritten  as 
(23')  P(r,t)    =  P^-plV^ 


-cA 

2 

(a-c)r^        2r'- 

.         ^2 

-(a  -  c)r 

'     o/  •          a2   2 
2(a  -  c)    r     . 
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In  (25)  the  argument  of  g  and  of  a  is  x,  given  by  (2J4.),   iVhen 
changes  which  occur  during  the  time  —  ,  which  is  that  required 
for  a  sound  wave  to  travel  one  bubble  r^'-dius,  can  be  neglected, 
{2k)    yields 

(26)  x  =  t-f  .ii^.t-f    . 

For  large  values  of  c,    where  only  the  first  terra  in  (17) 
need  be  retained  for  g(t),  (25)  becomes 


(27)   P(r,t)  -Pq-  dL^T"^ 


2  -2 
Po  -  PL  ■!■ 


-(a  a)  ^  i(a  a)  ]  _^_   (a  a)  _^ [(a  a)  ] 
Lr(l--)   2r-'(c-a;      2r  ■    ^r>.-c) 


.  2-.-    .  2..2n 
■(a  a)   _^  (a  a) 

"  ^       2r^ 


The  second  expression  in  (27)  has  the  same  form  as  th:>.t  given 
by  the  incompressible  theory  ( c  =  00  )  but  here  the  fvuiction 
a(t)  is  different  from  that  of  the  incompressible  theory. 
Also  the  argument  of  a  r.nd  a  is  x  instead  of  t,  shot^ring  that 
the  pressure  pulse  spreads  uith   finite  velocity. 

The  first  term  in  either  bracket  above  is  important  when 
a  is  changing  rapidly,  v:hich  occurs  near  bubble  minima,  and 
thus  is  mainly  responsible  for  the  pressure  in  the  primary 
shock  wave  and  the  subsequent  pressure  peaks.   The  factor 

(1  — )  in  the  denominator  of  this  term  in  the  first  form  of 
(27)  may  increase  this  term  somewhat  over  the  value  it  would 
have  in  the  incompressible  theory,  where  this  factor  is  unity. 
The  last  term  in  the  above  brackets  is  due  to  the  motion  of 
the  water  and  is  called  the  "af tcrf low"  pressure.   It  is  most 
important  between  pressure  peaks,  or  between  bubble  minima, 

I4..   The  Dir'^ction  Fi'.^ld;  n.-rgilta 

To  determine  the  qualitative  behavior  of  the  solution 
v(a)  of  (21)  we  will  nov;  examine  the  direction  field  in  the 
v,a  plane  determined  by  the  right  side  of  (21).   Wc  first  note 
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that  this  fiels  has  just  cv;o  singular  points--a  saddle  point 
at  a  =  0,  V  =   _-,  and  a  spiral  point  at  a  -  1,  v  =  0. .  We  next 
observe  that  v'(a)  bocoiiies  infinite  on  the  lines  a  =  0,  v  =  0 
and  V  =  c  and  it  becordes  zero  on  the  line  h(v,a)  =  0,   It  also 
changes  sign  on  crossing  those  lines.   This  information 
sxiXfices  to  determine  the  topological  character  of  the 
direction  field.   In  order  to  present  a  more  precise  sketch  of 
this  field,  vre  vjill  first  cxaninc  the  curve  h(v,a)  =  0. 
The  equation  h(v,a)  =  0  mas'-  bo  rewritten  as 

(28)   a^Y  =  2(3Y-l)(v-:r-^)[(v-c)^-(2+3c^)(v-c)  -iic-2c^]"^   . 

The  denominator  has  only  one  real  root,  and  thus  a  becomes 
infinite  for  only  one  value  of  v,  provided  that 
-(|+c^)^+  (2c+c^)^>0.  This  is  satisfied  if 

c  >7t(/T-4)  =  .103,  which  we  assume  to  be  the  case.   The  root 

^  '                                            8 
is  then  easily  seen  to  lie  in  the  interval  3c<v<3c(l  + ^)  . 

27c'' 

Furthermore  a,  given  by  (28)  has  only  one  maxiinuxi  when  y  <   -^  , 

2  h  ^ 

VJith  this  information,  assvming  that  c  >  .103  -■^'-^  "^-^Y^Ij  ^ho 

direction  field  has  the  topological  behavior  indicated  in 

figT-tre  1.   VJe  see  that  all  solutions  with  initial  velocity  v=0 

and  0<a  <1  increase  to  a  luaximum,  then  decrease  to  a  ininimvim, 

and  perform  damped  oscillations  approaching  a  =  1,  y  =0.   This 

is  exactly  the  behavior  observed  in  practice.  V/hon  c  is 

infinite  (19)  becomes  the  equation  of  the  sirdPle  incompressible 

theory,  and  has  only  a  vortex  point  at  v  =  0,  a=l.   The 

direction  field  for  this  case  is  shovm  in  figure  2.   All 

solutions  are  periodic,  undamped  or;cillations. 

V7e  have  determined  a(t)  for  a  particular  set  of  Initial 

data  by  numerical  integration  of  (16).   The  result  of  this 

calc\;lation  is  shown  in  fi,r;urc  3  along  vrith  a  set  of  iioa3iu?ed 

values  of  the  bubble  radius  obtained  experimentally.   Ihosc 

experimental  points  v;oro  taI:on  from  [1],  p.  271,  fig,  8.1,  and 

apply  to  a  .55  lb  tctryl  charge  fired  at  300  feet  below  the 

svirfacc  of  water. 
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To  carry  out  the  integration  wo  must  icnov/  A  ,  P^,  K,  y>  P 

end  C,  or  an  equivalent  sot  of  data,  see  cqs.  (7)  and  (l8). 

From  the  experimental  data  uo   estimated  that  the  ultimate  bubble 

radius  in  our  case  would  be  6  inches  which  wo  thci-'ofore  chose 

as  the  unit  of  length  L.  .  The  first  maximum  bubble  radius  ;-;as 

measured  to  be  18  inches,  Vjhich  becomes  3  in  tci^ms  of  our  unit 

of  length.   Since  the  iiieasuj?ed  data  was  poor  initially,  we 

started  the  numerical  integration  from  the  first  maximum  using 

initial  velocity  0  and  initial  radius  3>  integrating  both 

forwards  and  backwards,  in -time  from  this  point. 

•J 
The  values   y  ~   1»25>    P   =   1   g/cm-^,    C   =   l,if85  n./sGC  wore 

vised.     These   yield, 

c  =  U6.5 

T  =  I4..33  milliseconds   , 

Prom  the  figure  it  seems  that  the  theory  agrees  satisfactorily 
with  the  observations.   It  is  in  particular  worthwhile  to  note 
tliat  the  calculated  first  oscillation  period  is  27.9  luilli- 
seconds  as  compared  vrith  the  reported  experimental  value  of  28 
milliseconds.   Some  of  the  discrepancy  in  the  radius  time  curve 
is  undoubtedly  duo  to  error  in  estimating  the  ultimate  radius 
from  the  experimental  curve,   Tliis  not  only  affects  the  unit 
of  length  but  also  the  initial  radius. 

Using  the  computeJ.  values  of  a(t),  the  functions  f(^)  and 
£^  {K)   were  computed  froii  (17).   Graphs  of  these  functions  arc 
shovm  in  Figures  I4.  and  'S   res-cctiv&ly.   These  functions  can  be 
used  to  determine  the  pressure  F(r,t)  which  is  given  in  (23). 
In  fact,  for  largo  distances  r,  the  cxcoss-prcssiirc  vs,  time 
curve  is  essentially  a  multiple  of  the  f'  curve.   For  very  small 
distances,  however,  the  other  terms  in  oq.  (23)  may  become 
important,  especially  between  the  bubble  minima.  Yet  at  r  =  6, 
i,o,  at  a  distance  of  two  maximum  bubblt.  radii,  the  correction 
due  to  the  neglected  tei  •in  is  at  most  .35  percent  vrtiich  occurs 
at  ^  =  1.0. 
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5 •   Small  Oscillati on s 

If  the  initial  radius  is  nearly  equal  to  1,  waich  is 
the  equilibriuin  radius,  we  raay  write 

(29)  a(t)  =  1  +b(t) 

If  we  assiome  that  b(t)  is  snail,  insert  (29)  into  (16)  and 
retain  only  linear  terms  in  b,  we  obtain 

(30)  b  +  ll  b  +  3Yb  =  0 

Thus 

(31)  b(t)=b(0)e^  cos  (J3Y  -  (^)^  t) 

The  solution  (31)  yields  exponentially  danped  sinusoidal 
oscillations.   During  one  period  b  decreases  by  the  factor 


-jrj^ 


e      ,  if  (4^)   is  small  compared  to  3y. 

The  energy  in  the  bubble  consists  of  internal  energy 
and  potential  energy  and  is  given  by 

1         . 

(32)  E  =  ^  A3[P^-K.il^  (.^^jl  a3)-Y]  =  p^(^)r  a3(l+.^i_:^)   . 

jl.et  us  call  the  excess  energy  the  difference  between  the 
actual  energy  and  that  corresponding  to  the  equilibrivun 
radius  a  =  1.   Then  fo:?  small  oscillations,  using  (2;)  in 
(32),  the  excess  energy  in  given  by 

1 

(33)  E  =  ?A-^)'^  '-^f   o2(t) 

From  (31)  Jind  (33)  we  neo  that  the  excess  energy  decreases 

2T-J37 
by  the  factor  o    ^        dv.ring  each  period. 

In  Table  I,  arc  prosontod  some  calcylr-tions  based  on  the 
small  oscillation  theory  just  described.   The  dimcnsionlcss 

sound  speed  c  =  wF^  p  was  computed  for  throe  depths  lO'J,  200, 
and  lOCO  feet  respectively.   The  sound  speed  C  =  Ik'^S 
meters/sec  and  density  p  =  1  gm/cm'^  were  used. 
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Table   I 
Calculations   Based   on  Small    Oscillation  Theory 

Y  =   1.25 


Dimenslonless 
sound  speed 
c 


Decrease  in 
amplitude 
■oer  period 


Number  of 
periods  for 
amplitude 
to  be  halved 


Period 
j  correction 
I  factor 


6 .   Conclusion 

The  present  theory  of  a  spherical  bubble  talies  account 
of  the  compressibility  of  the  water.   Consequently  it  yields 
damped  radial  oscillations  of  the  bubble  and  finite  propagation 
speed  for  the  press\ire  pulses  in  the  water,  in  agreement  v;ith 
observation.   The  theory  is  essentially  the  same  as  that  of 
Zoller  [3],  from  which  it  differs  in  minor  details.   The 
feature  which  makes  the  theory  simple  analytically  is  the  use 
of  the  wave  equation  for  the  water  motion. 

Zoller  computed  the  pressure  in  the  water  and  compared  it 
with  more  exact  calculations  based  on  namorical  integration  of 
the  non  linear  fluid  dynamic  equations  in  the  water.   These 
calculations  were  also  based  on   the  assumption  of  uniform 
tDrcssurc  through  the  bubble.   The  agreement  was  excellent 
near  the  bubble,  which  was  the  only  region  considered.  He 
also  compared  the  radius -time  curves  obtained  by   these  tv;o 
different  calculations  and  thoy  also  agreed  wdl,  but  the 
time  over  vrhich  they  were  extended  was  not  very  long. 
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Our  calculations  constitute  another  check  on  this 
theory.   V/e  have  computed  the  radius- time  curve  for  nearly 
four  oscillations,  and  have  found  it  to  check  very  well  with 
an  experimental  curve  corresponding  to  the  same  initial 
conditions.   Thus  we  may  conclude  that  in  water  this  theory 
is  adequate  to  describe  the  bubble  oscillations  and  the 
pressure  in  the  water  near  the  bubble. 

It  is  particularly  important  to  point  out  that  in 
these  calculations  all  the  energy  is  accounted  j?or.   Since 
the  radius-time  ciorve  checks  with  experiment,  we  conclude 
that  all  the  energy  in  the  experimental  situation  is 
similarly  accounted  for.   Thus  the  question  of  "missing 
energy",  which  was  raised  in  recent  years  seems  unnecessary. 

It  should  also  be  pointed  out  that  certain 
compressibility  effects  v/ere  taken  into  account  by  Herring  [Lj.] 
in  an  attempt  to  estitiate  energy  loss  from  the  bubble.   More 
recently  Trilling  [5]  also  considered  compressibility  of  the 
water  in  showing  that  wave  motion  v;ithin  the  bubble  did  not 
modify  the  bubble  motion  daring  the  first  half-period  from 
that  computed  on  the  basis  of  uniform  pressure  throughout 
the  bubble . 

3y  examining  the  period  for  the  first  oscillation,  we 
rind  that  it  does  not  differ  appreciably  from  that  rriven  by 
the  incompressible  theory  since  c  is  large.   Therefore 
calculations  of  the  psriod  based  on  the  income  res sible  theory, 
such  as  those  of  Priodrrvan  [6],  can  still  be 
employed  for  the  first  period.   The  subsequent  periods  will 
duffer  because  of  energy  losses  due  to  comnrossibility. 
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Appendix;   Derivation  of  the  V/ave  Equation 

Vie   consider  the  non-viscous  motion  of  a  corapressible 
fluid  in  v;hich  pressure  is  a  function  of  density  alone. 
The  equations  of  motion  for  such  a  fluid  are 

(Al)  P^  +  V-  (pu)  =  0 

(A2)  u^+  (u  .V)u+  p"^^F  =  -  ^^ 

(A3)  p  =  p(p)   . 

In  these  equations  p  is  the  density,  p  the  pressure,  u  the 
velocity  vector,  and  Xi-  the  potential  of  external  forces  per 
ujiit  mass.   V/e  further  assum.e  that  the  motion  is  irrotational, 
which  implies  the  existence  of  a  velocity  potential  (j)  such  that 

{Ak)  u  =  V({)   . 

Prom  (A2),  (AL.)  we  have  the  Bernoulli  equation,  in  which 
F(t)  is  arbitrary: 


(A5)         ^(  V<I>)^  +  f   f^  +  i^  +  A  -  F(t)   . 


Now  differentiating  (A5)  with  respect  to  t  yields 
(a6)      V't>  •  V^t  "^  P""'"Pt  "^  '''tt  ^  -^t  "  ^t 

With  the  use  of  (A3),  the  derivatives  of  p  appearing  in  (Al) 
may  be  eliminated  in  terms  of  derivatives  of  p.   These  in  turn 
may  be  eliminated  by  means  of  (A2)  and  (A6).   We  have,  in  fact 

(A7)      p-^p^  =  p-^p^Pt  =  Pp[Ft  -  -^t  -  *tt  -  ^^  •  ^^t  ^ 


If  pressure  is  a  function  of  density  alone,  the  flow  is 
called  barotropic.   An  isentropic  flow  is  alv/ays  barotropic. 
Thus  if  the  entropy  is  initially  constant,  and  if  viscosity 
and  heat-conduction  are  absent,  the  flow  is  isentropic  and 
thus  barotropic,  in  the  p.bscncc.  of  shocks. 
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(a8)    p"^Vp  =  d"-^d.,VF  =  D^[-V-fL-  u,  -  (u.^)u]   . 
Using  (A7),  (A8)  in  (Al)  yields 
(A9)   K7^<t>  -  Pp<t>tt 

Equation  (A9)  is  the  exact  non-linear  equation  for  the 
velocity  potential,  in  v;hich  p   is  to  be  clininated  by  means 

Ox  (A3)  and  (a5).   When  the  right  side  is  negligible  and  p 

,  P 

constant  we  see  that  ^   satisfies  the  linear  wave  equation 

_  1  /2        , 
with  sound  sneed  c  =  p  '"'^  .   If  p  <Pj.4.  is  also  neglip:iblc,  9 

satisfies  Laplace's  equation. 

In  the  present  problem  JL =  0,  F.  can  be  set  equal  to 

zero,  and  the  terms  remaining  on  the  right  arc  derivatives 

,2 
of  ^ I     .       If  we  assume  the  speed  (^  ^)    is  small  compared  to 

c 
sound  speed  c,  these  terms  can  be  neglected,  yielding  the 

wave  equation.   VJe  also  assume  that  p   is  constant  which 

seems  valid  for  the  range  of  pressures  we  consider  when  the 

medium  is  water. 
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